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Abstract 



Coannihilation effects in neutralino relic density calculations are examined 
for a full range of supersymmetry parameters including large tan (5 and large 
Aq for stau, chargino, stop and sbottom coannihilation with the neutralino. 
Supergravity models possessing grand unification with universal soft breaking 
(mSUGRA), models with nonuniversal soft breaking in the Higgs and third 
generation sparticles, and D-brane models with nonuniversal gaugino masses 
were analysed. Unlike low tan/3 where too is generally small, stau coanni- 
hilation corridors with high tan/3 are highly sensitive to Aq, and large Aq 
allows too to become as large as ITeV. Nonuniversal soft breaking models at 
high tan (5 also allow the opening of a new annihilation channel through the s- 
channel Z pole with acceptable relic density, allowing a new wide band in the 
too — TO1/2 plane with m\ji ~ 400 GeV and too rising to 1 TeV. The D-brane 
models considered possess stau coannihilations regions similar to mSUGRA, 
as well as small regions of chargino coannihilation. Neutralino-proton cross 
sections are analysed for all models and it is found that future detectors for 
halo wimps will be able to scan essentially the full parameter space with 
m i/2 < 1 TeV except for a region with \i < where accidental cancellations 
occur when 5 ~ tan (3 ~ 30. Analytic explanations of much of the above 
phenomena are given. The above analyses include current LEP bounds on 
the Higgs mass, large tan/3 NLO correction to the b — > 57 decay, and large 
tan /3 SUSY corrections to the b and r masses. 
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I. INTRODUCTION 



Supersymmetry (SUSY) models with R-parity invariance offers a leading candidate for 
the dark matter observed in the universe at large and locally in the Milky Way galaxy Thus 
in models of this type, the lightest neutralino, Xv * s generally the lightest supersymmetric 
particle (LSP) and hence is absolutely stable. The observed dark matter is then the relic 
neutralinos left over after the Big Bang. SUSY models are predictive since they can calculate 
simultaneously the amount of relic density expected, Milky Way detection signals and the 
production cross sections at accelerators. Thus cosmological, astronomical and accelerator 
constraints simultaneously constrain the SUSY parameter space. 

The procedure for calculating the relic density of neutralinos is well known. (For a 
review, see M.) Over the past decade, a number of refinements of the analysis needed to 
get accurate answers have been included. Thus the treatment of threshold and s-channel 
resonances in the annihilation cross section in the early universe was discussed in |fj |J. 
In calculating the relic density, in most of the parameter space it is sufficient to use the 
two body neutralino annihilation cross section, since the effects of heavier particles are 
suppressed by the Boltzman factor. However, in special situations, the next to lightest 
supersymmetric particle (NLSP) may become nearly degenerate with the LSP, and the 
coupled annihilation channels of LSP-LSP, NLSP-LSP, NLSP-NLSP must be simultaneously 
considered . This phenomena of coannihilation was considered within the framework of the 
minimal supersymmetric model (MSSM) in || for the situation when the lightest chargino, 
Xi, becomes nearly degenerate with the Xi- in this paper we consider supergravity (SUGRA) 
models with grand unification of the gauge coupling constants and soft SUSY breaking at 
the GUT scale Mq = 2 x 10 16 GeV. The low energy properties of the theory is thus obtained 
by running the renormalization group equations (RGE) from Mq to the electroweak scale, 
where SUSY breaking at M G triggers SU{2) L x U{1)y breaking at the electroweak scale. 
Aside from being in accord with the LEP data implying grand unification, radiative breaking 
of SU {2)lxU (l)y greatly enhances the predictiveness of the theory as it reduces the number 
of SUSY parameters, and one finds generally for such models that the Xi is mostly Bino, and 
the xt is mostly Wino (when the soft breaking masses are less than 1 TeV). Then the xi 
and Xi do not become nearly degenerate and this type of coannihilation does not generally 
take place. (An exception for a class of D-brane models is discussed below.) 

More recently it was pointed out that in SUGRA models, the sleptons (particularly the 
lightest stau, f\) can become nearly degenerate with the Xi leading to a new type of coanni- 
hilation, and this was explored for low and intermediate values of tan/? = (H 2 ) / {Hi) 0,0. 
(Here H 2 ,i give rise to (u,d) quark masses). In this paper we examine this effect for the 
full range of tan/3, i.e. 2 < tan/3 < 50 M. Further, in addition to f\ — Xi coannihilation, 
we find the possibility of light stop t\ — Xi coannihilation as well as light sbottom b\ — x\ 
coannihilation. 

We consider these effects for a range of models based on grand unfication of the gauge 
coupling constants: (1) Minimal supergavity model (mSUGRA) M where there is universal 
soft breaking masses at Mq, (2) Nonuniversal supergravity models flIU |, where the first two 



generation soft breaking masses are kept universal at Mq (to suppress flavor changing neutral 
currents) as well as the gaugino masses, but the Higgs and third generation soft breaking 
masses are allowed to be nonuniversal, and (3) a D-brane model JTT|] based on Type IIB 
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strings |L2| where soft breaking masses at Mq of the SU (2) l doublet squark, slepton and 



Higgs are kept degenerate but distinct from the SU(2) L singlets, and similarly the gaugino 
SU (2) L doublet soft breaking mass is distinct from the SU (2) L singlets. While the first 
two models have been characterized as SUGRA models, they can also be realized in string 
models, and in fact any string model based on grand unification at a high Mq scale with 
the Standard Model gauge group holding below Mq, phenomeno logically can be treated as 
a SUGRA model with an appropriate amount of nonuniversal soft breaking. Thus the three 
models sample the possibilities of universal soft breaking, nonuniversal soft breaking in the 
squark, slepton and Higgs sector but universal gaugino masses, and finally nonuniversality 
in the gaugino sector as well. 

Since the Milky Way is perhaps 90% dark matter, it is a "laboratory" for studying the 
properties of dark matter (DM). Possible signals for DM include annihilation in the halo 
of the Galaxy, annihilation in the center of the Sun or Earth, and direct detection from 
scattering by terrestial nuclear targets. Of these, the last is most promising, and we restrict 
our discussion here to this. (For recent discussions of the other possibilities see [p~3| P~5T ) . 
In general, the neutralino-nucleus scattering has a spin independent and spin dependent 
part. However, for heavy nuclei (aside from exceptional situations discussed below), the 
spin independent scattering dominates, giving rise to approximately equal scattering by 
neutrons (n) and protons (p) in the nucleus. It is thus possible to extract from any data 
the neutralino-proton cross section, u^i-v ( su bject to astronomical uncertainties about the 
Milky Way). Current detectors (DAMA, CDMS, UKDMC) are sensitive to cross sections 

a^o_ p ~ 1 x l(T 6 pb (1) 

with a one to two order of magnitude improvement possible in the near future. Future 
detectors (e.g. GENIUS, Cryoarray) plan to be sensitive down to 

a -o_ p ^(10- 9 -10- 10 )pb (2) 

and thus it is of interest to see what parts of the SUSY parameter space can be examined 
by such detectors. 

In order to obtain accurate calculations of both the relic density and o"^o_ p for large 
tan f3, it is necessary to include a number of corrections in the analysis, and we list these 
here: (i) In relating the theory at Mq to phenomena at the electroweak scale, the two loop 
gauge and one loop Yukawa renormalization group equations (RGE) are used, iterating to 
get a consistent SUSY spectrum, (ii) QCD RGE corrections are further included below 
the SUSY breaking scale for contributions involving light quarks, (iii) A careful analysis of 
the light Higgs mass is necessary (including two loop and pole mass corrections) as the 
current LEP limits impact sensitively on the relic density analysis, (iv) L-R mixing terms 
are included in the sfermion (mass) 2 matrices since they produce important effects for large 
tan p in the third generation, (v) One loop corrections are included to m& and m T which 
are again important for large tan/3, (vi) The experimental bounds on the b — > 57 decay 
put significant constraints on the SUSY parameter space and theoretical calculations here 
include the leading order (LO) and NLO corrections for large tan f3. Note that we have not 
in the following imposed b — r (or t — b — r) Yukawa unification or proton decay constraints as 
these depend sensitively on unknown post-GUT physics. For example, such constraints do 
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not naturally occur in the string models where SU(5) (or 50(10)) gauge symmetry is broken 
by Wilson lines at Mq (even though grand unification of the gauge coupling constants at 
Mq for such string models is still required). 

In carrying out the above calculations we have included the latest LEP bound on the 



light Higgs mass, > 114 GeV [16[] and the recent NLO corrections for large tan/? for 
the b — > 37 decay [|i~7l , |i~8f . (We have checked numerically that [17j and [18|] give identical 
results.) Since there are still some remaining errors in the theoretical calculation of rrih 
as well as uncertainty in the top quark mass, m t = (175 ± 5) GeV, we will conservatively 
assume here that rrih > 110 GeV. In the MSSM, the constraint on rrih is tan/? dependent as 
Ah production (A is the CP odd Higgs) when rriA = Mz is possible for tan f3 ~ 9 and this 
can be confused with Zh production. However, in SUGRA models, radiative electroweak 
breaking eliminates this part of the parameter space, and so it is correct to impose the LEP 
bound on m h for the full range of tan/3. LEP also gives the bound |19| m-± > 102 GeV, 



and the Tevatron bound for the gluino (g) is ]20|] m- g >270 GeV (for squark (q) and gluino 
masses approximately equal). We assume an allowed 2a range from the CLEO data for the 
b — > sj branching ratio |2T[ : 

(3) 



1.8 x 10~ 4 < B(B -> X sl ) < 4.5 x 10" 
In calculating the relic density we will assume the bounds 



0.02 < Sl^k 1 



where VL^o 



< 0.25 



(4) 



Px°/pc p x ° is the relic density of the Xi, Pc = 3Hq/8itGn (H is the present 
Hubble constant and Gn is the Newton constant) and Hq = /il00kms _1 Mpc _1 . The lower 
bound on Q x »h 2 is somewhat lower than more conventional estimates (Q x oh 2 > 0.1), and 
allows us to consider the possibility that not all the DM are neutralinos. (In the following, 
we will mention when results are sensitive to this lower bound.) Accurate determinations 
by the MAP and Planck satellites of the dark matter relic density will clearly strengthen 
the theoretical predictions, and already, analyses using combined data from the CMB, large 
scale structure, and supernovae data suggests that the correct value of the relic density lies 
in a relatively narrow band in the center of the region of Eq. (|4]) | 2"2fl . We will here, however, 
use the conservative range given in Eq. (|4]). 

Supersymmetry theory allows one to calculate the x?-q uar k cross section and we follow 
the analysis of p| to convert this to Xi ~ V scattering. For this one needs the tcN <j term, 



o-kN = 7^(m u + m d )(p\uu + dd\p), 



(5) 



— (m u + m ( i)(p\ss\p) and the quark mass ratio r = 2m s /(m u + m^). We use 
30 MeV H, and r = 24.4 ±1.5 |25|. Recent analyses, based on new 7r — N 

wn 



O"0 — O'ttN 

here a = 

scattering data gives cr nN 



We 



65 MeV @. Older n — N data gave a nN ^ 45 MeV _ 
will use in most of the analysis below the larger number. If the smaller number is used, 
it would have the overall effect in most of the parameter space of reducing <J x o^ p by about 
a factor of 3. However, in the special situation for fi <0 where there is a cancellation of 
matrix elements, the choice of a^^ produces a more subtle effect, and we will exhibit there 
results from both values. Some of the results described below have been mentioned earlier 



in conference talks 28 . 
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There have been in the recent past a number of calculations of cr^o_ p in the literature 



24,29-34 



However, none include simultaneously all the corrections discussed above, the 
large tan (3 range and the recent LEP bounds. We find general numerical agreement with 
other calculations for regions of parameter space where the omitted corrections are small. 



II. MSUGRA MODEL 

We consider in this section the mSUGRA model. This model is the most predictive of 
the SUGRA models as it depends on only four additional parameters and one sign. These 
may be chosen as follows: m , the universal scalar soft breaking mass at Mq; m^, the 
universal gaugino mass at M G (or alternately one may use m^o or m g as these quantities 
approximately scale with m^, i.e. m^o ~ 0Ami/ 2 and m g ~ 2.8m!/ 2 ); A , the universal 
cubic soft breaking coupling at Mq] tan/? = (H 2 ) / (Hi); and n/\n\, the sign of the Higgs 
mixing parameter fi in the superpotential term W = fiHiH2. We restrict the parameter 
space as follows: 

m ,mi/ 2 < 1 TeV (6) 
|A)/mi /2 |<4 (7) 
2 < tan/5 < 50 (8) 

The relic density is calculated by solving the Boltzman equation in the early universe. If 
the excited SUSY states lie significantly above the x?, then this takes the standard form [|TJ: 

$ = ~3Hn - (a ann v)(n 2 - n 2 eq ) (9) 

where n is the neutralino number density, n eq is its equilibrium value, (...) means thermal 
average, a ann is the annihilation cross section Xi + Xx ~ > f (f — fi na l state), and v is the 
relative velocity. In this case, the effects of the heavier particles are suppressed by the 
Boltzman factor. If however, one or more particles, Xi, have masses mi near the neutralino 
(Xi = Xi, m^o = mi), coannihilation effects can occur. In this case Eq. @ holds for the 
neutralino number density with, however, n eq and cr ann given by 



n, 



cq 



J2n ieq (10) 



J2 r i r 3 a W r i = -T^ ( U ) 



i,j n eq 



In the nonrelativistic limit where freezeout occurs, the equilibrium number densities take 
the form P,|35| 



n ieq = 9l (^) 3 / 2 exp(-m t /T)[l + ^] (12) 

Here crjj is the cross section for Xi + Xj — > / {a\\ is the neutralino annihilation cross section) 
and gi are the degeneracy factors. In our analysis, we have included all final states that are 
energetically allowed. 
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One sees from Eqs. (|TT|JT^ ) that the effects of coannihilation depend simultaneously on 
the near degeneracy of any heavier particle with the neutralino, and on the size of the 
annihilation cross sections. A large cross section can then compensate partially for the 
Boltzman factor. In SUGRA models two effects occur which enhance the coannihilation for 
the sleptons. First, since the neutralino is a Majorana particle, its annihilation cross section 
is p-wave suppressed. Consequently, 

Cij, <Jn = O[10<7n]; i, j = sleptons (13) 

and the slepton cross sections are enhanced relative to the neutralino cross sections an. Sec- 
ond, a numerical accident allows the selectrons to be nearly degenerate with the neutralino 
over a region of parameter space of increasing mo, mi/2- One can see this analytically for 
low and intermediate tan j3 for mSUGRA where the RGE can be solved analytically |36| . 
At the electroweak scale one has 

K B = m l + (6/5) hm\ /2 - sin 2 9 W M^ cos(2/5) (14) 



m-o = {a x /a G )mx/2 (15) 

where = [1 - (1 + (3 t t)- 2 }/ ft, t = \n(M G /M z ) 2 and ft is the 17(1) (3 function. Numerically 
this gives for e.g. tan (3 = 5 

m\ R = ml + 0.15m 2 /2 + (37GeV) 2 

m|o = 0.16m 2 /2 (16) 

Thus for mo = 0, the cr becomes degenerate with the Xi a t m i/2 — 370 GeV, i.e. coanni- 
hilation effects begin at = (350 — 400) GeV. For larger m^, the near degeneracy is 
maintained by increasing m , and there is a corridor in the m — mi/2 plane allowing for an 
adequate relic density. (This is in accord with the numerical calculations of [|3C|j.) 

For larger tan/5, the situation is more complicated, as the the L-R mixing in the f\ 
(mass 2 ) matrix, generally makes the f\ the lightest slepton (lighter than the e# or fin), and 
hence the fi dominates the coannihilation phenomena. The RGE must of course now be 
solved numerically, and there is strong sensitivity to A . 
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FIG. 1. Allowed corridors in mSUGRA in the mo — m^M plane satisfying the relic density 
constraint for fi > 0, tan f3 = 40 for (from bottom to top) A = m 1 / 2 , 2 m l / 2 , 
terminate at low mj^ due to the b —* sj constraint. 



Ami The corridors 



The new effects are exhibited in Fig. 1 where the allowed regions in the m — tTixii plane 
satisfying Eq. @ for A =m 1 / 2 , 2m 1 / 2 and 4mi/ 2 are shown for /i > and tan/? = 40. One 
sees several features of the large tan /3 region there. First, the allowed corridors are sensitive 
to the value of Aq, the allowed values of itlq increasing with Aq. Further, the corridors end 
at a minimum value of mi/2 due to the b — > 57 constraint. (We note that these lower bounds 
on rrii/2 are sensitive to the NLO corrections to the b — > 37 branching ratio mentioned in 
Sec. 1 above.) In fact, there is only the coannihilation region left in the parameter space, 
since the b — > 57 constraint at large tan (3 has eliminated the non coannihilation part of the 
parameter space. (Actually, there can be small islands of allowed parameter space left in 
the non coannihilation region (mi/2 ~ 350 GeV) for large tan/3 due to a cancellation that 
can occur between the Standard Model amplitude and the chargino amplitude.) While for 
low tan/3, m is generally small ( m o ~ 200 GeV), we see that at large tan/3, mo can 
become quite large reaching up to 1 TeV for large Aq. Finally we mention that while for 
fixed A the allowed corridors are relatively narrow, as one allows A to vary from small to 
large values, the allowed region due to coannihilation covers most of the m — mi/2 plane. 

We turn next to consider the Xi~P cross section. Large cross sections can arise for large 
tan/3 and small m^o, and in [33|, it was shown that Eq. ([I]) could be satisfied in mSUGRA 



only also for small Q^oh 



1. e. 



for 



tan/3 ~ 25, Vt^h 2 ~ 0.1 



Xl 



m*o ~ 90 GeV 



(17) 



Small cross sections generally arise for small tan/?, large m and large mi/ 2 . It is of interest 
to see what the minimum cross section predicted by theory is in order to know how much 
of the SUSY parameter space will be examined by proposed detectors. In order to see the 
effects of coannihilation, we first examine the domain mi/2 < 350 GeV {m^o < 140 GeV) 
where coannihilation does not take place. Fig. 2 exhibits the minimum cross section for 
tan/? = 6, fi > 0, satisfying all the accelerator and relic density constraints. One sees that 
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9 pb for tan (3 > 6, and so the entire parameter space in this region would 
be accessible to detectors planned with sensitivity of Eq. (0). 



CTvO_„ ~ 1 x 10 
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FIG. 2. 



(<7^o_ p ) m i n in the noncoannilation region as a function of m^o for mSUGRA for fj, > 0, 



tan /3 = 6, obtained by varying over the range of mo and Aq. 



At higher m.1/2, coannihilation occurs, allowing mo and m^j to become larger, and 
hence the cross section should become smaller. We examine first the case /i > 0. The 
A dependence at large tan/5 shown in Fig. 1 then implies that the cross section should 
also decrease with increasing Aq, since, as seen in Fig. 1, mo then increases. This is seen 
explicitly in Fig. 3. 
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FIG. 3. o"^o_ p as a function of mi/ 2 for mSUGRA fj, > 0, tan/3 = 40 and Aq = 2m 1 / 2 (upper 
curve), 4m 1/ / 2 (lower curve). 

Thus it is possible that as mi/ 2 gets large, a consequence of coannihilation is that large 
tan j3 and large A can give smaller cross sections than low tan f3. This is illustrated in Fig. 
4 where cr^o_ p is plotted as a function of m^ for A = 4m 1 / 2 , tan/3 = 40 (upper curve) and 
tan/3 = 3 (lower curve). One sees that for rrii/ 2 > 600 GeV, the tan/3 = 40 curve actually 
drops below the tan (3 = 3 curve, the A dependence (producing large m ) compensating for 
the tan f3 dependence. One has however, the bound 



<Jv°-n ~ 1 x 10 



-10 



pb, for /i > 0, m-i/2 < 1 TeV 



(18) 



S 



Again, most of this parameter space will be within the reach of planned future detectors. 
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FIG. 4. cr^o_ p as a function of m 1 / 2 for mSUGRA, fj, > for Aq = 4m 1 / 2) tan/? = 40 (upper 
left curve) , tan (3 = 3 (lower left curve) . 

We turn next to the case of fi < 0. The situation here is more complicated. As pointed 
out for low and intermediate tan/? in ||31|| , an accidental cancellation can occur between 



the Higgs w-quark and the Higgs d-quark scattering amplitudes which can greatly reduce 
c^o-p and we investigate here whether this cancellation continues to occur in the high tan/? 
region. Since the cancellation is somewhat subtle, we briefly describe how it comes about. 
In general, Xii is a mixture of Bino(-B), Wino(VK) and the two Higgsinos (Hi, H 2 ): 



Xi = N n B + N l2 W + NizHi + N U H 2 (19) 

where Nu are the amplitudes. In much of the parameter space, the t-channel Higgs exchanges 
(h, H) dominate cr^o_ p . The d and u quark Higgs amplitudes are 



~ 2M W \ cos (3 m\ cos (3 m 2 H ) 
A u = 92 m u ( cos a F h sin a F H \ 

~ 2M W \sm (3 ml sin Pm 2 H J [ ' 

where a is the Higgs mixing angle and 

F h = (N 12 - N n tan 9 W ) (N u cos a + N 13 sin a) (22) 
F H = (Ni 2 - Nu tan 9 W ) (N u sin a - N 13 cos a) (23) 

Since the s-quark contribution to the scattering is quite large, the (i-quark amplitude A d 
will generally be quite large. However A d will be suppressed if the amplitudes AT 13 , obey 
the equation 



tan a + cot a 

r u ~ -a 13 (24) 



In general tana is negative and small (tana = — cot (3 ~ —0.1), and Eq. ([24]) can happen 
if A14/A13 is positive, which is the case for ji < 0. Once A d is sufficiently suppressed, it can 
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cancel the remaining (smaller) w-quark contribution coming from A u , leading to a nearly 
zero value of cr^o_ p . For fixed tan/3, total cancellation occurs at a fixed value of m 1( / 2 , though 
the effects of cancellation has a width. What happens is shown in Fig. 5 where cr^?-p * s 
plotted in the large TO1/2 region for tan/3 = 6 (short dash), tan/3 = 8 (dotted), tan/3 = 10 
(solid), tan/3 = 20 (dot-dash) and tan/3 = 25 (dashed). One sees that the cross section 
dips sharply when tan/3 is increased from 6 to 8, and and goes through a minimum for 
tan/3 = 8 at mu2 = 810 GeV. For tan/3 = 10, the minimum recedes to mu2 = 725 GeV, 
and then begins to advance for higher tan/3 i.e. rising to m\/2 = 830 GeV for tan/3 = 20, 
and rai/2 = 950 GeV for tan/3 = 25. Thus if we restrict mu2 to be below 1 TeV, the cross 
section will fall below the sensitivity of the planned future detectors for mi/2 ~ 450 GeV for 
a restricted region of tan (3 i. e. 

a^o_ p < 1 x 10" 10 pb for 450 GeV < m 1/2 < 1 TeV; 5 ~ tan/3 ~ 30 (25) 
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FIG. 5. cr-o_ p for mSUGRA for /i< 0, A = 1500 GeV, for tan/3 = 6 (short dash), tan/? = 8 
(dotted), tan/3 = 10 (solid), tan/3 = 20 (dot-dash), tan/3 = 25 (dashed). Note that the tan/3 = 6 
curve terminates at low due to the Higgs mass constraint, and the other curves terminate at 
low mi/2 due to the b — ► 57 constraint. 

At the minimum, the cross sections can become quite small, e. g. < 1 x 10~ 13 pb, without 
major fine tuning of parameters, corresponding to an almost total cancellation. (Actually, 
as pointed out in |[38||, the true minimum of the cross section would then arise from the 



spin dependent part of the scattering, which though very small, does not possess the same 
cancelation phenomena.) Further, the widths of the minima for fixed tan/3 are fairly broad. 
While in this domain the proposed detectors obeying Eq. (fj) will not be able to observe 
Milky Way wimps, mSUGRA would imply that the squarks and gluino would lie above 1 
TeV, but at masses that would still be accessible to the LHC (i.e. < 2.5 TeV). Also note 
that this phenomena occurs only for \i < and for a restricted range of tan /3, so cross checks 
of the theory would still be available. 

The minima occurring in Fig. 5 arose from cancellations between the Higgs w-quark and 
Higgs <3>quark amplitudes, and as such are sensitive to the quark content of the proton. 
As discussed in Sec. 1, cr^/v is 110 1 yet we h determined, and Fig. 5 used cr nN = 65 MeV. 
In Fig. 6, we compare this choice(solid curve) for tan /3 — 10 with the parameters used in 
27l , |3lH (dashed curve), where a n N = 45 MeV. One sees that the minimum at mi/ 2 = 725 
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GeV is shifted to m\ii = 600 GeV, with analagous shifts occurring for the other values of 
tan (3. Thus the extreme cancellations that occur in the coannihilation region for /x < are 
quite sensitive to the properties of the proton. 
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FIG. 6. o-o_ p for mSUGRA for p < 0, tan/3 = 10, A = 1500 GeV for a nN = 65 MeV (solid 
curve) and the parameters of [26,36] (dashed curve) where ovjv = 45 MeV was used. 



III. SUGRA MODELS WITH NONUNIVERSAL SOFT BREAKING 

We consider here models with nonuniversal soft breaking masses. In order to suppress 
flavor changing neutral currents, we maintain the universal soft breaking slepton and squark 
masses, m , at Mq in the first two generations. However, experiment is compatible with 
nonuniversal Higgs and third generation soft breaking at Mq, and these can arise from 
nonuniversal couplings in the supergravity Kahler potential. We parameterize this situation 
as follows: 

m Hx = m o( x + <h); m H 2 = m o( 1 + $0; 

TO ff J =mg(l + <f 3 ); m t 2 R = mg(l + 5 4 ); m T 2 R = mg(l + S 5 ); 

m 2 R = m 2 (l + 6 6 )\ m 2 = m 2 (l + 8 7 ). (26) 

where = squarks, II = (^ t >Tl) sleptons etc. The Si measure the deviations from 

universality. (If one where to impose SU(5) or 5*0(10) symmetry then 5 3 = 5 4 = S 5 = S w , 
and Sq = S-j = S5.) In the following we limit the Si to obey 

-KSi<+l (27) 

The lower bound is necessary to prevent tachyons at Ma, and the upper bound represents 
a naturalness condition. We maintain gauge coupling constant unification at Mq (which is 
in accord with LEP data). We also assume here gaugino mass unification. (Deviations from 
this can arise from nonuniversal couplings in the supergravity gauge function f ab , and an 
example of this is treated in the next section.) 

While the nonuniversal models contain a number of additional parameters, one can obtain 
an understanding of the new effects these imply from the following analytical considerations. 
In the decomposition of Eq. (|19D , the lightest neutralino is mostly Bino, i.e. iV n ~ 0.8. The 
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Higgs mixing parameter /i to a large extent controls this mix, and as fi 2 decreases (increases) 
the higgsino content increases (decreases). In order to see qualitatively what happens we 
examine the low and intermediate tan/5 region where the RGE can be solved analytically. 
One finds for /i 2 the result (see e.g. |39|| ): 



t 2 



1-3A) 



t 2 - 1 

1 + A> 



S 2 + 



2 

Si 

t 2 



1 

t 2 * 



I -D 



(6 3 + S 4 ) 



m 2 , + universal parts + loop corrections. 



(28) 



where t = tan/5, and Dq = 1 — (m i /200GeVsin / 9) 2 = 0.25. One sees that the universal 
contribution to the wig term in fi 2 is small, making fi 2 quite sensitive to nonuniversal con- 
tributions. Also since Do is small, the Higgs and squark nonuniversalities enter coherrently, 
approximately in the combination S% + 5^ — 5 2 . (The RGE can also be solved analytically for 
the SO(10) model where tan/3 > 40, and numerically otherwise, with results qualitatively 
similar to Eq. (P8f).) In addition we note here the mass of the CP odd Higgs boson, itla is 



in 



A 



t 2 +l 

l + Do 



3(1 



2 

S 2 +61 



D r 



(S3 + S4 



m\ + universal parts + loop corrections. 



(29) 



The effects of the nonuniversalities can be seen to fall into the following categories: 

(i) The spin independent neutralino-proton cross section depends upon the interference 
between the gaugino and higgsino parts of the neutralino, the larger the amount of inter- 
ference the larger the cross section. Thus CTyO_ p can be significantly increased relative to 
mSUGRA by decreasing fi 2 (and hence increasing the higgsino content) by chosing <5 3 , <5 4 , 
Si < 0, 5 2 > 0, and be decreased (though by not so large amount) with the opposite sign 
choice. 

(ii) In the early universe annihilation cross section, the s-channel Z° pole amplitude 
depends on the higgsino contribution of the neutralino (the Bino and Wino parts giving zero 
coupling). Thus if /i 2 is lowered by chosing 83, 5 4, Si < 0, 5 2 > one gets an increased amount 
of annihilation, and a decreased amount with the opposite choice of signs. Thus deviations 
from universality can significantly effect the satisfaction of the relic density constraints of 
Eq. (§• 

(iii) The third generation squark and slepton nonuniversalities shift the mass spectrum 
of the particles. Thus the positions of the corridors of coannihilation, and their widths can 
be significantly modified (also modifying the cross section cr^o_ p .) 



Effects of type (i) were early observed in |39| , and further discussed in where is was 
seen that cr^-p could be increased by a factor of 10 or more with the sign choice 53,4,1 < 0, 
5 2 > 0. We here update this result to take into account of the latest LEP bound rrih > 114 
GeV. Fig. 7 shows the maximum value of c^-p f° r the nonuniversal models for fi > with 
the above choice of signs for 5j for tan/3 = 7 (lower curve) and tan/3 = 12 (upper curve). 
One sees that detectors obeying Eq. (JTJ) are sampling the parameter space with tan/3 > 7 
(compared with the mSUGRA result of tan/3 > 25). 
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FIG. 7. Maximum value of a^o_„ as a function of m,~o for the nonuniversal model with a > 0, 
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To illustrate further some of the effects of nonuniversality, we consider first the case 
where there are only Higgs mass nonuniversal soft breaking, i.e. only 5 12 are non zero. 
In general for the universal mSUGRA case, both /i 2 and rriA are determined in terms of 



the other SUSY parameters. From Eqs. (|28|) and (|29"D, we see that these quantities now 
depend on the additional parameters 5\ and 62 , and thus may be varied keeping the other 
(mSUGRA) parameters m , mi/ 2 , A and tan/5 fixed. In fl4"0|| , it was thus assumed that fi 2 
and rriA could be chosen arbitrarily. Actually, however, both fi 2 and m^ are constrained by 
the bounds of Eq. fl2"7|), and varying /1 2 and m^ arbitrarily can lead to points in parameter 
space with unreasonably large values of 6^2 (or even tachyonic values with <5i j2 < — 1). Thus 
here we will examine the actual parameter space subject to the conditions of Eq. ([27]) (as 
well as the usual relic density and accelerator bounds). As an example of an effect of type 
(ii), we consider in Fig.8 the case in the coannihilation region where A Q = m\/2, tan/3 = 40, 
and chose 82 = 1. Since fi 2 is insensitive to 5i we set Si = 0. In this domain, the allowed stau 
coannihilation region corresponds to the lowest band in the mo — mi/2 plane of Fig. 1, with 
mo running from about 200 GeV to 400 GeV and the corresponding corridor is only slightly 
effected by 62 7^ 0. As shown in Fig. 8, the effect of the nonuniversal soft breaking now 
allows the opening up of a new wide band of allowed region at much higher values of m . 
Thus normally such large values of mo would give rise to too little annihilation to satisfy 
Eq. (^). However, by lowering /i 2 and hence raising the higgsino content of the neutralino, 
the Z° channel annihilation is enhanced allowing for an acceptable relic density. In Fig. 8, 
the low side of the band corresponds to fi^/i 2 =0.25, while the upper end to Q^oti 2 =0.02. 
If one were to raise the relic density lower bound to Q^oh 2 =0.1, the width of the band 
would be reduced by about 30% to 40%. (The small dip'at m 1/2 = (400 - 450) GeV (i.e. 
m^o =(160 - 180) GeV) arises due to the opening of the t — t channel in the annihilation cross 
section.) Thus we see that nonuniversal soft breaking allows for much larger values of mo 
than mSUGRA and hence a heavier mass spectrum, and this could be an experimental signal 
for nonuniversalities. (Universal soft breaking only allows for the coannihilation corridor in 
Fig. 8 at much lower m .) In addition, the larger m leads to a smaller o"^,o_ p . 
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FIG. 8. Effect of a nonuniversal Higgs soft breaking mass enhancing the Z° s-channel pole 
contribution in the early universe annihilation, for the case of 82 =1, tan/3 = 40, Aq = mi/2, 
\i > 0. The lower band is the usual fi coannihilation region. The upper band is an additional 
region satisfying the relic density constraint arising from increased annihilation via the Z° pole 
due to the decrease in /j, 2 increasing the higgsino content of the neutralino. 

We consider a second example where both effects of type (ii) and (iii) occur simultane- 
ously, by examining third generation squark and slepton nonuniversal masses obtained when 
(5io(= S3 = 64 = 65) takes on the value 5iq = -0.7. (All other 5i are chosen zero.) We again 
assume A = m 1( / 2 , tan/3 = 40 and fi > 0. The effect of 5 W is shown in Fig. 9. The lowest 
band is the usual mSUGRA f\ coannihilation band (shown here for reference only). The 
upper band corresponds to two phenomena: The lower half is the actual f\ coannihilation 
region which has been shifted upwards due to the fact that m 2 , for the f% has been replaced 
by mo(l + 5io), and since 5\q is negative, one needs a larger mo to achieve the coannihilation 
corridor. The upper part corresponds to the opening of the Z° channel as in Fig. 8 since 
5 W < reduces /j 2 increasing the higgsino content of the neutralino. For mi/2 ~ 500 GeV, 
the two bands overlap. However at higher m^, they are separated. Thus in the high m\M 
region, the low side of the middle band (where f\ coannihilation is occuring) corresponds 
to tt^oh 2 =0.02, while the upper side corresponds to ^o/i 2 = 0.25. In the region between 
the two top bands, fl^oh 2 >0.25 (and hence excluded). Finally, since 5 3 and 5 4 are negative 
(reducing fj, 2 and also making the total m 2 , contribution to /i 2 negative), for sufficiently large 
tuq, the Z° channel opens reducing ^^oh 2 to 0.25 on the lower side of the top band, that 
band terminating when f^o/2 2 = 0.02 on the top. 
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FIG. 9. Allowed regions in the rao-m^ 2 plane for the case tan (3 = 40, A$ = m^, \i > 0. The 
bottom curve is the mSUGRA f\ coannihilation band of Fig. 1 (shown for reference). The middle 
band is the actual f\ coannihilation band when <5io = —0.7. The top band is an additional allowed 
region due to the enhancement of the Z° s-channel annihilation arising from the nonuniversality 
lowering the value of /j, 2 and hence raising the higgsino content of the neutralino. For mj/j ~ 500 
GeV, the two bands overlap. 

We note in both Figs. 8 and 9, that the effect of nonuniversalities can greatly enhance 
the allowed region in the m — m\/2 plane for large m^, particularly by increasing m , and 
these changes occur without unreasonable amounts of nonuniversality. 

We consider next the neutralino-proton cross sections arising in these nonuniversal soft 
breaking cases. One might have expected that the new regions in the m — plane allowed 
by the relic density constraint would correspond to smaller values of cr^o_ p since they are 
bands of much larger mo- However this is compensated by the fact that the nonuniversal 
effects considered above also reduce /i 2 and effects of type (i) compensate for the largeness 
of mo and in fact give a net increase to the cross section. This is illustrated for the example 
of Fig. 8 in Fig. 10 where cr^o_ p is plotted for 5 2 =1, tan/3 = 40, A = m\/2- The solid 
line corresponds to the cross section for the f\ coannihilation corridor (the lower band of 
Fig. 8), while the long (short) dashed curves correspond to the cross sections along the top 
(bottom) of the Z° pole band. We see that the Z° band cross sections, where mo is very 
large, are quite large, and this region of parameters should soon be accesible to e.g. to the 
next phase of the CDMS experiment in the Soudan mine. The cross section in the lower mo 
domain corresponding to the f\ coannihilation region are much smaller, and to explore this 
domain would require e.g. the GENIUS or Cryoarray detectors. 

In Fig. 11, we examine the expected neutralino-proton cross section for the case of 
<5io = —0.7, tan/3 = 40, A = mi/2, > of Fig. 9. Here the large dashed curve is the 
cross section running along the top of the upper band of Fig. 9 (the top of the Z° channel 
band) and the short dashed curve is the cross section along the bottom of the middle band 
(the bottom of the f\ coannihilation corridor). The solid curve is shown for comparison 
and would be the cross section in the f\ coannillation corridor for universal soft breaking 
(mSUGRA). We see again that the Z° corridor with very large mo would be accesible to 
the CDMS Soudan experiment. The actual f% coannihilation corridor cross section lies 
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considerably higher than the universal case, but would probably require the GENIUS or 
Cryoarray detectors to explore. 
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FIG. 10. cr^o_ p as a function of m^ 2 for 5 2 = 1 (other Si = 0), for tan/3 = 40, Aq = m 1 / 2 , 
[i > 0. The long (short) dashed curves are the cross sections along the upper (lower) sides of the 
Z° bands of Fig. 8. The solid curve is cr^o_ p along the fi coannihilation corridor of Fig. 8 . 
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FIG. 11. a x±-v as a f unc ti° n °f m i/2 f° r £10 = -0-7 (other <^ = 0), for tan/3 = 40, Aq = mi/2, 
> 0. The long dashed curve is the cross section along the top of the upper band in Fig. 9 
where Z° annihilation is occurring. (The reduction at m.1/2 ~ 500 GeV is due to the b — ► S7 
constraint.) The short dashed curve is the cross section along the bottom of the middle band 
where f\ coannihilation effects occur for this nonuniversal model. The solid curve is included for 
comparison and is the cross section in the f\ coannihilation corridor for the case of universal soft 
breaking (mSUGRA). 

In the above discussion, we have chosen the signs of the b% parameters to decrease /i 2 
to illustrate the dramatic effects such 5i can produce without any unusually large amount 
of nonuniversality. If instead one were to chose the opposite signs, e.g. 5 2 < 0, or 5 W > 0, 
the fj 2 would have been increased. For this region of parameter space the Z° annihilation 
channel would now be absent (as in mSUGRA). The usual f\ coannihilation corridors would 
of course still be present, but for the 5 W > case at lower values of m . The cross sections 
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for n > would then be somewhat reduced, but still accessible to future detectors. 

The cross sections discussed above were for /x > 0. We now consider the \x < case. In the 
mSUGRA model, it was seen that for large m\/2 with 7711/2 < 1 TeV and a restricted range of 
tan/3, a special cancellation could occur which reduced cr^-p to be below 10~ 12 pb (see Fig. 
5), well below what future detectors could see. This cancellation involved the Higgs w-quark 
and Higgs (i-quark amplitudes, when the neutralino components obey Eq. (plj) . However, 
nonuniversal choices such as those of Figs. (8-11), involve now two region of parameter space 
consistent with the relic density bounds: the coannihilation region, and the region above 
it dominated by the s-channel Z° pole. Eq. Q23j ) can in fact hold for the coannihilation 
corridors, and the minima seen in the mSUGRA models occur here, at essentially the same 
points as in Fig. 5. However, for the Z° bands, /i 2 is quite reduced, and Nxa is increased 
violating Eq. (|24|). Thus no such minima occur for mo in the Z° bands of Figs. 8 and 9. 



IV. D-BRANE MODELS 



Recent studies of Type IIB orientifolds have shown that it is possible to construct a 
number of string models with interesting phenomenological properties |12| . The existance of 



open string sectors in Type IIB strings implies the presence of Dp-branes, and in the simplest 
cases, one can consider the full D = 10 space compactified on a six torus T 6 . Models of this 
type contain either 9-branes and 5;-branes, % = 1, 2, 3, or 3-branes and 7j-branes. Associated 
with each set of n coincident branes is a gauge group U(n). There are a number of ways 
in which one can embed the Standard Model in such systems. We consider here the model 
where SU(3)c x U(l)y is associated with one set of 5-branes, 5i, and SU(2)l with a second 
intersecting set 52 [0. Strings starting on 5i and ending on 52 then carry the joint quantum 
numbers of the two sets of branes, the massless modes then being the quark and lepton left 
doublets, and the two Higgs doublets. Strings starting and ending on 5i then carry the 
quantum numbers of that brane, the massless modes being then the quark and lepton right 
singlets. 

Supersymmetry breaking can be treated phenomenologically by assuming that the F 
components of the dilaton S and the moduli Tj grow VEVs f!2]]4"I|] : 



F s = 2V3 < ReS > sm6 b e iaa m 3/2 (30) 
F T > = 2V3 < ReTi > cos^e.e^mg/a (31) 

where 9 b , 0j are Goldstino angles (O 2 + O 2 . + 2 = 1). In the following we assume 9 3 = 0. 
T-duality determines the Kahler potential |T2|j41[| , and using Eqs. (30,31) one generates the 



soft breaking terms at Mq for the gauginos, 

fh l = m 3 = — A = v / 3cos6 , 60ie _Ml m 3 /2 (32) 
m 2 = ^cos^l -0 2 ) 1/2 m 3/2 (33) 

and for the squarks, sleptons and Higgs, 



12 - (1 - 3/2 sin 2 6 b )m 2 3/2 for q L , l L , H u H 2 (34) 
m 2 = (1 - 3 sin 2 6» 6 )m 2 /2 for u R , d R , e R . (35) 
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Originally this model was considered to examine the possible effects of the new CP violating 
phases in Eqs. (32,33) and in the \i and B parameters |TT| , |42[| . However, it was subsequently 



seen that unless tan /3 is small (i.e. tan (3 ~ 5) an unreasonable fine tuning of the GUT scale 
parameters is needed if the experimental constraints on the electron and neutron electric 
dipole moments are to be obeyed ||42|| . The LEP data has now eliminated a considerable 



amount of the low tan (3 region, and since we are interested here in effects that arise for large 
tan f3, we will in the following set all the phases to zero. The model then depends on three 
parameters: 9b, Oi, and m 3 / 2 , with 9b < 0.615 and < ©i < 1, and contains a unique type 
of nonuniversal soft breaking. In terms of the notation of Eq. (26) one has 

5 l = 5 2 = 6 3 = 5 7 = -(3/2) sin 2 9 h (36) 
5^ = 5$ = 5§ = —3 sin 2 9b (37) 

where we have equated mo to m 3 / 2 . Note that these nonuniversalities hold for all three 
generations, not just the third generation. In addition the gaugino masses at Mq are no 
longer universal (except when 0! = 1/a/2), and there are two "effective " gaugino masses at 
Mq which we label by m\/2 = fh\ = rh 3 , and my 2 = m^- The A$ parameter in this model 
is then fixed by A = -m^. 

We consider first the question of what parts of the parameter space might be accessible 
to current detectors which obey Eq. (1). From Eqs. (34,35) we see that the sfermion masses 
decrease with increasing 9b, and so we expect the largest cross sections to arise for small 9b, 
and of course for large tan f3. However, the b — > 37 constraint eliminates the higher values 
of 9b- The limiting boundary obeying Eq. (1) is shown in Fig. 12, where cr^o_ p is plotted 
for tan (3 = 20 and 9b = 0.2. For this case, the neutralino mass would be quite small but is 
still consistent with the LEP bounds^]. (The gap in the middle is due to over annihilation 
in the early universe through the s-channel Z pole.) The cross section rises, of course with 
higher tan (3, and so current detectors are sampling regions where 

tan/3~ 25 (38) 

It is easy to see that slepton coannihilation effects exist in this model similar to those 
of the SUGRA models. Since both the cr and Xi both evolve from the GUT scale using 
the U(l) gaugino, Eq.(16) still holds with m 2 , replaced rriy 2 (l — 3 sin 2 9 b ). Thus for fixed 9b, 
one can adjust 777,3/2 so that the Xi an d are nearly degenerate, allowing for coannihilation 
to take place. In addition, however, there can be regions of parameter space where there is 
chargino-neutralino coannihilation. This arises because of the nonuniversal gaugino masses 
at Mq. Thus from Eqs. (32,33), when 0i gets large, m 2 becomes smaller than mi, and the 



LEP determines the bound on m^o by measuring lower bounds on the chargino mass and lower 
bounds on m^o + m^o |^] . This combined with lower bounds on the Higgs and slepton masses and 
the value of the Z width, limits the SUSY parameter space giving rise to the m^o bound. Thus 
the bounds depend on the SUSY model used, and LEP quotes 37 GeV for the MSSM, and 48 GeV 
for mSUGRA. We have checked that for the D-brane model being considered, the LEP data still 
implies m^o > 37 GeV. 
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Xi becomes mainly Wino, and can become degenerate with the Xi ■ Coannihilation becomes 
possible, however, only for a very small range in ®i near 0.9. 
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FIG. 12. a^o __ for D-brane model for u > 0, 6^ = 0.2 and tan/3 = 20. The gap in the curve 
is due to excessive early universe annihilation through the s-channel Z pole. 



As in SUGRA models, a cancellation of matrix elements can occur for fi < 0, allowing 
the cross sections to fall below the sensitivities of planned future detectors. This is exhibited 
in Fig. 13, where cr^o_ p is plotted for tan/3 = 6 (solid curve), 12 (dot-dash curve), and 20 
(dashed curve). (The tan/3 = 6 curve terminates at low m^o due to the constraint, while 
the higher tan/3 curves terminate at low m^o due to the b — > S7 constraint. The upper 
bound on m^o, corresponding to rrig < 1 TeV, arises from the Qh 2 constraint.) One sees 
that the cross section goes through a minimum at tan/3 = 12, though the expected rise at 
higher m^o does not appear since the parameter space terminates before this sets in. 
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FIG. 13. a Xi-p ^ or ^ e D-brane model for fi < and tan/3 = 6 (solid), tan/3 = 12 (dot-dash), 
and tan/3 = 15 (dashed). 

V. STOP AND SBOTTOM COANNIHILATION. 



It has been known for some time that the nearness of the t-quark Landau pole can 



significantly lower the i\ mass ||44|| , and this raises the possibility of i\ — x coannihilation 
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occurring [f|5[ . Working against this possibility is the positive contributions to m~ from the 
ST/ (3) gaugino contribution, and of course the m 2 contribution. Nonuniversal soft breaking, 
however, opens additional possibilities of this type, as it can allow the coefficient of wig to 
become negative, reducing further the t\ mass. To examine this in more detail, we consider 
the low and intermediate tan (3 regime where the RGE can be solved analytically. For the 
Ir one finds 

l- A), 



m 



tR 



A) + ( 
(1- 



-(5 2 + 5 3 

o 

)D ^A m 1/2 
r 



H o — o 4 



m n 



1 - D )D A 2 + em\ /2 



h 



3/3) ^>I/2 + m t + 3 sin Q wMz cos 2(3 (39) 

The fi form factors, i = 1,2, 3, are numerically approximately 40, 100, 700 respectively, and 
e, and H3/F given in |4l|, are approximately -2 and +2 respectively. One sees that for Aq 
large and negative, the Landau pole A Q terms do give a significant negative contribution 
and m~^ would be minimized for (5 2 3 > and 64 < 0. The universal Wq term, scaled by 
Do = 0.25, is already quite small, so it does not require too much nonuniversality to achieve 
a degeneracy with the Xv ( in fact this can be done even for the universal case if Ao and rrii/2 
are chosen sufficiently large. However, the degeneracy occurs in that case for parameters 
in the TeV domain.) Since early universe annihilation cross sections involving the stop are 
strong interaction phenomena compared with the electroweak neutralino cross sections, the 
coannihilation effect can be quite significant. However, it is necessary to impose the usual 
constraints, particularly the b — > 37 bounds, and also that the stop must lie below the stau. 
Thus a satisfactory coannihilation region occurs only when mo is large (e.g. greater than 
about 800 GeV). The neutralino-proton cross section in this domain generally lies above 
10~ 10 pb, and thus should be accessible to future detectors. 

Stop coannihilation phenomena are more subtle in the D-brane model, since the nonuni- 
versalities are related i.e. 5 2 = 63 = —(3/2) sin 2 9bmy 2 = 5^/2. The wiq term for this case 
is 



[D - (l + 2Do)sm 2 e b ]m 



3/2 



(40) 



which is negative for sin 9 b ~ 1/6. (There are also modifications in the gaugino terms 
to acount for their nonuniversal nature in the D-brane model.) Thus stop coannihilation 
effects can occur also in this model. Here, however, the parameter space for coannihilation 
is narrowed and generally occurs when m 3 / 2 > 1 TeV. 

6-squark coannihilation phenomena can also occur in these models. For low and inter- 
mediate tan (3, where one may neglect the 6-quark Yukawa coupling (relative to the t-quark 
coupling), the &r mass is generally large (due to the SU(3) gluino contribution). However, 
the bi mass can get small for nonuniversal soft breaking. Thus rag is given by 
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and the light 6-squark can be achieved by chosing 5 2 , 5 4 > and £3 < so that the coefficent 
of the Wq term is negative. In this case the light 6-squark is mostly bi, and since the t 
and b squarks form an SU(2) doublet, the light stop lies just above the the light sbottom, 
and there can be regions where both contribute simultaneously to coannihilation. However, 
these coannihilation regions again occur for mo > 1 TeV. 



VI. CONCLUSIONS 

We have considered in this paper coannihilation effects in dark matter neutalino-proton 
detection cross sections which arise for large tan/? and large values of Aq for models which 
have grand unification of the gauge coupling constants at the GUT scale Mq. Coannihi- 
lation can occur if the light stau (fi), light chargino (xf), light stop (ii) or light sbottom 
(61) become nearly degenerate with the neutralino. In order to examine these effects for 
a wide range of models, we have considered SUGRA models with universal soft breaking 
(mSUGRA), nonuniversal soft breaking in the Higgs and third generation squarks and slep- 
tons, and finally a D-brane model based on type IIB orientifolds which possess nonuniversal 
gaugino masses. In all these we have imposed the cosmological neutralino relic density 
constraints as well as accelerator constraints on the SUSY parameter space (including in 
particular the large tan /3 NLO corrections to the b — > 37 decay flT7| , |IB| as well as the large 
tan/5 SUSY corrections to m& and m T ). Radiative breaking of SU{2) x U{\) at the elec- 
troweak weak scale also gives strong constraints on the theory. 

mSUGRA illustrates some of the effects arising for large tan /3. Thus for low tan /3, the 
regions in the mo — plane allowed by the relic density constraint Eq. (fj) is relatively 
insensitive to A . However, for large tan/3, the b — > 37 constraint allows only the T\ — %\ 
coannihilation corridors to remain (except for some small islands when ~ 350 GeV) 
and as A increases, these corridors have increasing m . Thus m can rise to 1 TeV for large 
mi/2- (See e.g. Fig. 1.) This lowers the value of cr^o_ p and even though u^o_ p increases 
rapidly with tan /3, the value of 0"^o_ p at large tan (3 and large Aq can be smaller than cr^o_ p at 
tan/3 = 3, as illustrated in Fig.4. For /1 > 0, however, one generally finds for the parameter 
space of Eqs. (6-8) that cr^o_ p ~ 1CT 10 pb, and so should be accessible to future detectors. 
In the high tan/3 domain, the b — + 57 constraint produces a lower bound on m^, which 
decreases with increasing A$. Thus one finds that m^o ~ 120 GeV (m§ ~ 740 GeV) for 
A = 4mi/ 2 , tan/3 = 40, \i > 0. 

For /i < 0, a special cancellation can occur in mSUGRA allowing cr^o_ p to become much 
less than 10~ 10 pb, as has been previously noted for low and intermediate tan/3 [|31 1 . The 



analytic origin of this effect is discussed in Sec. 2. The effect, however, occurs only for a fixed 
range of tan /3. Thus the position of these minima begins for tan /3 = 7 at = 1000 GeV, 
moves to lower m^j until tan/3 increases to 10, and then increases to beyond mi/2 = 1000 
GeV at tan/3 ~ 30. (See Fig. 5.) We note that while the cancellation occurs at a fixed value 
of mi/2-, the effects spread over a wide range of m^j and in this region, halo dark matter 
would be inaccessible to current or planned detectors. However, if the recently reported 
anomaly in the muon gyromagnetic ratio |47j is indeed due to supersymmetry, then /i is 



positive and this case would be eliminated. 

The nonuniversal soft breaking models illustrate new phenomena at large tan/5 not 
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present in mSUGRA, and one can also understand these effects analytically. Thus the value 
of /x 2 , obtained from radiative breaking of SU (2) x £7(1), is a major element in controlling the 
value of &x°- P , and lowering ji 2 increases cr^o_ p . Thus for one sign of the nonuniversalities, 
cr%n-p can ^ e increased by a factor of 10 or more, allowing current detectors to sample the 
parameter space with tan/3 as low as = 7. (See Fig. 7.) For large tan/3 a new effect can 
occur. The value of fj, 2 also controls the gaugino/higgsino content of the Xi, and lowering \j? 
increases the amount of Higgsino in the Xi- This then increases the amount of early universe 
annihilation via an off shell s-channel Z boson, which opens an additional allowed region in 
the m — my 2 plane, considerably wider than the f\ coannihilation corridor and at much 
higher m . (See Fig. 8.) (Other simple choices of nonuniversal parameters can also raise the 
values of mo for the f\ coannihilation corridor, as illustrated in Fig.9.) Though mo is large 
for this Z-channel annihilation, reducing ct^o„ p , this effect is compensated by the fact that 
fi 2 is decreased (which raises a^o_ p ). Thus cr^o_ p remains relatively large (see Figs. 10, 11), 
and this region of parameter space should become experimentally accessible in the relatively 
near future (e.g. when CDMS moves to the Soudan mine). 

The D-brane model examined is quite constrained as all the squark and slepton nonuni- 
versalities are controlled by a single parameter 9b- Current detectors are sampling parts of 
the parameter space where tan (3 ~ 25. Because these models have also nonuniversal gaugino 
masses, the LEP lower bound on m^o is much smaller than in SUGRA models, i.e. the same 
as in the MSSM: m^o > 37 GeV. However, the gaugino nonuniversality still allows the usual 
type of fi — Xi coannihilation to occur. A unique feature of the gaugino nonuniversality is 
that 1712 can become smaller than ih\ allowing the Xi to become mostly Wino and a xf ~ Xi 
coannihilation domain to exist. However, the near degeneracy of xf and Xi occurs only for 
a very small region of parameter space (when Qi = 0.9). The cancellation phenomena in 
cr^o_ p for /i < also occurs for these D-brane models, the minimum value of m^o where the 
cancellation is complete occuring at m^o = 315 GeV for tan/3 = 12. 

Nonuniversal SUGRA models also allow for t\ — Xi coannihilation to occur. In this 
case, the light stop is mostly tf>. Then for Aq < 0, the Landau pole contributions to m 2 R 
are negative, and one can chose nonuniversal soft breaking so that the coefficient of m\ is 
also negative. The t\ then becomes nearly degenerate with the Xi when mo ~ 800 GeV, 
coannihilation setting in for that domain. Actually, t\ — Xi can also occur in mSUGRA, 
but to achieve the near degeneracy one needs A > Amy 2- Similarly, this coannihilation can 
occur for the D-brane model but for m 3 / 2 > 1 TeV. 

In a similar fashion, b\ — Xi coannihilation can occur for nonuniversal models when 
the b\ is mainly hi- Thus one may chose Aq < and nonuniversal parameters in m~ so 

that the coefficient of m\ is negative, making the b\ nearly degenerate with the Xv This 
coannihilation only sets in, however, for m > 1 TeV. 
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Note added: After completing this paper there appeared the paper "The CMSSM Parameter 
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Space at Large tan/?" by J. Ellis, T. Falk, G. Ganis, K.A. Olive and M. Srednicki, [hcp- 



ph/0 102098| . This paper considers large tan/3 effects for the mSUGRA model only, and is 
restricted to A = 0. 
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